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There are strong arguments favoring a four-quark interpretation of sub-GeV light scalar mesons 
and the diquark-antidiquark body-plan of the tetraquark seems to provide the most convincing 
picture. The building diquarks of these particles are assumed to be spin zero objects. In this paper 
we explore the possibility that radially excited aggregations of spin zero or spin one diquarks might 
exist and discuss the possibility of the y(2175) state observed by BaBar and confirmed by BES 
being one such state. 

'qq" ' Introduction. It has been shown recently how the diquark-antidiquark interpretation of the sub-GeV 

f^ ' scalar meson nonet made of /o(980), ao(980), k(800), (^(SOO) can lead to a remarkable description of the decay 

^^ ' properties of these particles [l[, adding a rather strong confidence that they are indeed tetraquark objects. 

CN ' In terms of diquarks q, the light scalar nonet can be interpreted as qq particles made of spin zero diquarks, so 

1— H , called "good" . A good diquark operator in the attractive anti-triplet color channel (greek letters) , antisymmetric 

, ^ ' in flavor (latin letters) can be written as [^] : 

^^ ■ q^a = eijk<^af3-(qi'^l5q'''^ (1) 

(-H Also "bad" , spin one, diquarks can be conceived. A bad diquark operator can be written as: 

J"; < = ea/3^(gt^7g'^ + q^Mqn (2) 

5-2-1, 

~5 Both represent positive parity, O'^ and 1+, states. Similarly one can construct 0^ and 1^ operators as: qcQ 

and Qc^jsq. The latter are identically zero in the 'single mode configuration', quarks that are unexcitcd with 

respect to one another. Lattice studies, see, e.g. 0], suggest that diquarks are preferably (energetically) formed 

in spin zero configurations. In fact the most solid tetraquark candidates are scalars made of good diquarks. 

The next step in building this new spectroscopy is to find states belonging to other nuiltiplets. Recently, 

the BaBar experiment has observed a new J^^ = 1 resonance, the y(2175), decaying into (f>fo [J|, later 

confirmed by BES Q . The fact that it has been first observed into a mode including a light scalar makes it a 

f*^ ' good candidate for belonging to a higher tetraquark multiplet. 

In this paper we make a simple Ansatz that properly reproduces the light scalar nonet and utilize it to make 

predictions on possible excitations, focusing on J^'~" — 1^^ states and searching for a match with the F(2175). 

QQ ' The Model. This is developed in the framework of a non-relativistic Hamiltonian including spin-spin 

^^ , interactions inside a single diquark, spin-spin interactions between quarks and antiquarks belonging to the two 

diquarks forming the hadron, spin-orbit and a purely orbital term: 



1^ 

o 



X 



H = 2m^ + i/(|') + H^f + H,, + H,, , (3) 



C^ where: 



^sV ^ ^Kg^g^ {Sq^ ■ Sq^ + Sq^ ' "S, J -f ^Kq^g^ Sq^ ■ Sq^ + "^Kq^q^ Sq^ ■ Sq^ 

H.=«,^^. (4) 

The parameters in these equations are fit to data: m^ is the mass of the [q^q.^] diquark, Kq is the spin-spin 
coupling between the quarks inside the diquarks, Kqi^j are the spin-spin couplings ranging outside the diquark 
shells, Aq is the diquark spin-orbit coupling, and B^ weights the contribution of the total orbital angular 
momentum of the qq system to its mass; the overall factors of two are just conventional notations. We focus 
on the case where only light u, d, s quarks are involved. The spin-spin interaction Hamiltonian has the form: 

pairs ■^ 



because the color-magnetic moments are inversely proportional to quark masses. In Eqs. ([4]) we incorporate the 
mass dependencies in the k^ constants. The Hamiltonian ([S]) describes contact interactions. For this reason we 
could expect that allowing a relative orbital angular momentum between the diquarks will decrease or switch-off 
the spin-spin interactions between quarks and antiquarks, namely H^'j,''' . In the following we shall consider both 
cases. 

The values of the couplings appearing in Eqs. Q were estimated in Ref. 6] from a fit to meson and baryon 
masses under the assumption that the spin-spin interactions are independent of whether the pair of quarks 
belong to a meson or a diquark. The estimates are summarized in Tab. [H 
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TABLE I: Estimate, in MeV, of the parameters in the Hamiltonian in Eq. |4] depending on the diquark type q = [giQa] 



Extending the same procedure to the S* = 1, i = 0, 1 meson states p{770), ai(1230), 02(1320), 6i(1229) we 
also infer the parameters related to the orbital angular momentum: A^ — 22.5 MeV, B^ — 505 MeV. 
To describe a qq quantum state we adopt the following non-relativistic notation: 



Pq, iq! Jqq/ — \s i^ Q, S i 9; <Jqq/ 



(6) 



where T oa a2 for a spin zero diquark and F ex ai for a spin one diquark. The action of a spin-spin interaction 
operator, e.g. Sg ■ Sq, on ^ is described as follows: 

(Ss ■ Sq)\s'^Tq, s'^Tq; S'qq) = ^ X! l*'^rCTj9, s'^crjTq; Sqq) (7) 

j 

As an example let us diagonalize the Hamiltonian in Eq. ([4]) between scalars made of diquarks, i.e. |0q, Oq; 0) 

(8) 



with a relative iqq — 0. With an obvious shorthand notation: 



|0q,0q;0) = -|s^cr2<?, s'^cr2g;0) := -a2<»<72- 



We can then compute: 



(0q,0q;0|5,-5,|0q,0q;0) = --x3, 



where we have used the fact that <jJ(J2 — ^^2<^j and <Jj<Jj =3x1. The final result is: 



m = 2r7lq — iKq. 



If q = [sq], then, using the values in Table [H we get: 

TO = 988 MeV, 



(9) 
(10) 
(11) 



reproducing the mass of TOqo and tti/q , considered as [qs] [qs] particles with the two diquarks in spin zero and in 
5"— wave [l|]. Repeating the same calculation with q — [ud] one gets, for the cr-meson mass: 



TO = 481 MeV 



(12) 



Higher Mass Tetraquark Spectrum. The next orbital excitation comes when Lqq = 1 and both good 
and bad diquarks are considered. Among these, also 1 multiplets are generated, which are the main interest 
of this paper. To estimate the masses, one needs to repeat the diagonalization with the basis: 



1) = |0q,0q;lj) 

|lq,Oq;lj) + |Oq,lq;lj) 



|2) 

|3) = |lq,lq;l. 



V2 



(13) 
(14) 
(15) 



Since both the good and the bad diquarks have positive parity, the state |2) has P = C = —I, provided that 



1. For the states |1) and |3), since Cqq(— 1) i^(— 1) 
1. 



1, c, 



qq 



-'m 
To perform the diagonahzation we adopt the shorthand notation described above: 

1 



-1 provided that S„ 



0, 2 and 
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e'^'=(T2^, 



2%/2 



0'20-fc, 



(16) 



Hence, it is rather straightforward to derive the mass term shift Arrigg due to the part of the Hamiltonian in 
Eq. ([4]) constraining only spin-spin interaction terms, H^^;'- 



Am, 



-3K,q 

Ka 



(17) 



^g "'9l?2 ''*9i?i +'*92«2''/ . 

Writing the latter matrix as diag(Ai, A2, A3), the four solutions for states having quantum numbers 1 are: 



,(1) 



TUy ('-'qi — U, Oq2 — U, Oqq — (J, ivqq — 1) ^ 2,Tflq^ + Ai + -Dq 

"^y (S'qi = l,5'q2 = 0, Sqij = 1,-^qq = 1) = StTIj, + 5 + A2 - 2A^ 



-Ben 






i, Oqj 
i, Oq2 



-L, *^qq — U, i^qq 



2m^ 



2(5 + A. 



2A, + B^ 



1, S'qq = 2, Lqq = 1) = 2mq + 2(5 + A3 - 6Aq + Bq 



(18) 



where 5 = m^(^s=i) ~ "^q(s=o)- Following Jaffe and Wilczeck i3], we will assume for q = [qq], 6 ~ 285 MeV 
whereas for q = [sq], 5 ~ 150 MeV. The numerical values for the coefficients of A^ and B^, call them a,b, are 
given in the folowing table. 





a(Sq, S^,L) 
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S^^l,Si,^0,L^l 


-2 




Oq = 1, Oq = 1, Oqq ^ Z, _L/ ^ 1 
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Oq ^ -L, Oq ^ -L, *-5qq --~ L^ ±j ^ L 


-2 




S'q =1,5^ = l,5'^ii =0,L = 1 


-2 





TABLE II: Eigenvalues of spin-orbit and angular momentum operators in l|4]l. All these combinations of diquark spins and orbital 
angular momenta allow a, J = L + S^^ = 1 state. 



In case q = [ss], only the last state in Eq. p6)) is allowed since only bad diquarks can be formed by Fermi- 
Dirac. One should therefore consider only the (Iq, Iq; lj|iJgg|lq, Iq; Ij) correction to the mass, from Eq. (fTT]) 
is equal to n^ — ingg. We therefore have: 



m^Y {Sqi = Ij'S'qa = 1,'S'qq = 0, Lqq = 1) = 2mq + 26 + {Ks - 2Kss) - 2Aq + Sq 
my (S'qj = l,5'q2 = 1,5'qq = 2, Lqq = 1) = 2mq + 26 + (Ks - 2Kss) ~ 6Aq + B^. 



(19) 

The numerical values for my'' masses can be found in Table Hill The fact that 1 tetraquark particles require an 
angular momentum barrier Lqq ^ between diquarks must suppress the diquark-antidiquark chromomagnetic 
interactions. Switching off spin-spin interactions between quarks and antiquarks leads to the second estimates 
in Table Uni 

Tetraquark Decay Modes. From these results, if we want to identify the F(2175) as a tetraquark state, 
we have to resort to cither a q = [qs] hypothesis, with spin-spin interactions between diquarks set to zero or 
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TABLE III: Mass values itIy in MeV for the 1 states as computed from Eqs. JTSj and II19I I. When apphcable, the first value 
includes spin-spin interactions between diquarks, the second one neglects them: H]?^ = 0. 



to a q = [ss] hypothesis: indeed we use the hypothesis that /o(980) in the decay products of y(2175) is itself 
a qq particle with q = [qs]. In order to test the match of the Y{2175) with these assignments, we study its 
possible decays under both hypotheses, q = [qs], [ss], and the correspondence with observations. With both 
assignments, the observed F(2175) -^ 0/o(98O) decay mechanism would be described by the following diagram: 



r(2175) 



S- 

s,q- 

s,q- 

s- 




q,s 



/o(980) 



As for other expected decays, a significant contribution is expected to come from (j)r] via the diagrams below: 



y(2175) 



s,q 



s,q 



(t^vW) 



Similar diagrams would also yield (fxj) and 7777 decays, but they are forbidden by charge conjugation and Bose 
statistics selection rules. 

We can also estimate the decay width of the F(2175) — > (t>ri channel. The decay proceeds through P-wave 
and the matrix element is given by: 



{^{p', 6(«)r7(q)|rb, e(^))) = gye'^'^P^p^qJ^^^e'^P 



(20) 



The quark exchange amplitude considered above has been studied first in Q were a rather good fit of the 
scalar meson decays to pseudoscalar, S -^ PP, was obtained associating to this amplitude a coupling strength 
A ~ 2.6 GeV. Discarding angular momentum barrier effects and following the definition given in [6[: 
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where here My = My, we get the following estimate for the Y partial width in (jjrj: 



r(y(2i75)^07?) = c: 



2 SttMI- 
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(21) 



(22) 



where ^ = g or | depending on the q — [qs] or q = [ss] respectively. We therefore estimate r(y (2175) -^ (prj) ^ 5 
MeV or ~ 20 MeV under the two hypotheses respectively. 

The most typical decay mode expected for tetraquarks is the baryon-antibaryon one. Stretching the color 
string between the diquark and the antidiquark a qq pair is formed, qOSSSSSS^iq^ZSSSSS^q, letting two baryons in 
the decay products. The favored decay diagram is 



r(2175) 



S- 

s,q- 

s,q- 

s- 



q,s 



S,A 



S,A 



where the topology is suggested by the fact that the diquark in the A baryon could only be of the [ud] type. 

Since the SS decay mode is phase-space forbidden for the y(2175), only the q = [sq] assignation would allow 
a dominant baryonic decay, y(2175) -^ AA. With the other assignment, q = [ss], the y(2175) -^ AA would be 
made possible by the annihilation of an ss pair by, e.g., an instanton interaction giving two pairs of light quarks 
in the final state: in diagrammatic terms 



The 6-fermion instanton interaction has the form Ci ex det(q^j^q-'j^) and its role in scalar meson dynamics has 
been recently underscored in [l|. 

Nonetheless it is known from fll] that the instanton coupling, fitted to explain light scalar meson decays like 
/o(980) — > TTTT, is about ten times smaller that the quark exchange one. Therefore, under the hypothesis that 
the baryonic mode would be instanton driven in the q = [ss] case, we would not expect it to be easily visible. 

The latter diagram allows also a number of possible decays of the y(2175) like pp, cttt, tttttt, and ■qiT^n'^ in 
the q — [sq] hypothesis and K^ n^ , K~ ir^ K'^ , and (J)tt°tt° in the q = [ss] one. 
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FIG. 1: Fit to the (fifo invariant mass distribution. 



Observed Y(2175) Decays. In order to test the compatibility of the F(2175) state with a tetraquark 
interpretation and to discriminate between the two possible diquark compositions, we have reanalyzed the 
published BaBar data for e+e^ -^ 4'fol 3, AA7 [lO|, and tfirij ll'|. These are initial state radiation processes, 
where J^^ = 1 states are produced together with the initial state photon. The invariant mass of the system 
produced with the photon is then expected to show a resonant behavior in correspondence to states. 

We perform simultaneous fits applying a consistent notation for the Breit-Wigner and several possible models 
for the non-resonant component A„r- The general notation for the expected cross section as a function of the 



invariant mass of the system under study is 

a{m) ex ^■^pg{m)\A{m)\ 
s 



(23) 



A{m) = e''A„r{m) + ^aoBiY^f) 



moTtotirno) 



m? -71111 + i^tot{m)mQ 



where 5 is the relative phase between the two components at the pole; ^pg{rn) = (p(77i)/p(mo))"^ is the final 
state dependent phase space factor: a^^ = 3, a</,/o = Oi/^x — 1; p{'m) is the momentum of the two particles in 
the final state when their c.o.m. energy is to; toq a-nd ctq are the pole mass and production cross section and 
are independent of the considered final state; B{Y — > /) is the Branching Fraction to the specific final state; 
Ttotl'Ti) is the comoving width, the sum over the considered final states plus a constant term to account for all 
other decays with thresholds far from the pole mass: 



rtot(TO) 



T,\l-Y.^{Y^f){l~e{m)) 



(24) 



e{m) = <^'ps{m)/FBw{m) 



where Fbw is the Blatt-Weisskopf factor [ij] and Fq is the bare width. Note that for masses below the threshold 
of a given final state the corresponding ^ is imaginary. 

TABLE IV: Fit results to the ^/o invariant mass distribution (Fit A), the ^/o and AA invariant mass distributions (Fit 
B), and the invariant mass distributions of all three modes (Fit C). The meaning of the symbols is explained in the text. 



Fit 


xVdof 


mo (MeV) 


Fo (MeV) 


R{AA,(t>fo) 


R{<t>ri,<t>fo) 


A 
B 
C 


25/36 
47/45 
85/103 


2167 ±11 
2158 ±11 
2153 ±9 


69 ±21 
66 ±20 
72 ±20 


N/A 
44 ±19 
6.6 ±3.5 


N/A 
N/A 
10 ±3 



Of the three considered modes, in two cases the y(2175) would decay to states where it is above threshold while 
for the other (AA) toq is below threshold. A fit to the discovery mode, i/i/o, with Anr = A x e^'^'(^~™*^™/o) x 
(1 - e-(^-"*-'"/o_)Vai) (see Fig.H]) returns the results listed as "Fit A" in Tab. |TVl 

Including the A A mode in a simultaneous fit, with common bare mass and width and letting the ratio between 
the branching fractions i?(AA, (/)/o) — Bp^j^/B^f,^ fioat, yields different results depending on the assumptions 
Anr in the AA mode (see Fig. [2]). In case we assume no non-resonant contribution, Anr = 0, we find that the 
K(2175) — > AA decay can explain the whole observed spectrum, with mass and width parameters consistent 
with the fit to the (f)fo mode. The results are reported as "Fit B" in Tab. IIVI 
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FIG. 2: Fit to the AA invariant mass distribution under the assumption of no non-resonant contribution (left), or 
assuming an exponential amplitude for it (right). 



There is an indication that the AA decay is favored, even by one order of magnitude. If instead we assume 
Anr{m) = A X e"*"'™, letting k float, we observe no significant decay into AA, but with huge uncertainties: 
R{AA,(t)fo) = 23 ±51. 

Finally, the (fu] mass distribution shows a significantly higher background than the other modes and small 
structure in a position which is lower and narrower than the one observed in the (pfo channel and marginally 
consistent with it. Fitting the three distributions simultaneously as shown in Tab. IIV| "Fit C" gives an overall 
good fit (see Fig.[3l). 




FIG. 3: (j>fo, AA, and (pr; invariant mass distributions with the simultaneous fit overlaid. 

We can then conclude that the presence of the F(2175) -^ AA is is suggested by the fact that its existence 
would explain the whole AA mass spectrum. Under this hypothesis the baryonic decay mode would be dominant 
{Bj^j^/B^fg = 44 ± 19), thus favouring q = [qs] for the F(2175). As for the (j>r] mass distribution, uncertainties 
are large, but the case of a relatively large decay amplitude into it is not disfavored. 

Both the calculation of the mass spectrum and the reanalysis of the experimental data tend to favor the 
assignment of the Y{2175) to a [s(7][sg] state with both diquarks in the S = 1 state and with one unit of relative 
orbital angular momentum. 

Conclusions . In this paper we have studied the consequences of allowing spin one diquarks to build 1 
(qq) orbitally excited tetraquark states potentially visible in processes with initial state radiation at BaBar 
and Belle. In particular we have focused on the y(2175) resonance recently discovered by BaBar. This particle 
could be the first tetraquark state showing the expected baryon-antibaryon decay. Indeed, reanalyzing BaBar 
data, we find that if we set to zero the non-resonant contributions, the AA decay mode is the prominent one, 
indicating a q = [qs] assignation for the F(2175). 
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